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Motivated by tlie disagreement found at liigli dissipation between simulation data for the heat 
flux transport coefficients and the expressions derived from the Boltzmann equation by the standard 
first Sonine approximation [Brey et al., Phys. Rev. E 70, 051301 (2004); J. Phys.: Condens. Matter 
17, S2489 (2005)], we implement in this paper a modified version of the first Sonine approximation 
in which the Maxwell-Boltzmann weight function is replaced by the homogeneous cooling state 
distribution. The structure of the transport coefficients is common in both approximations, the 
distinction appearing in the coefficient of the fourth cumulant 02- Comparison with computer simu- 
lations shows that the modified approximation significantly improves the estimates for the heat flux 
transport coefficients at strong dissipation. In addition, the slight discrepancies between simulation 
and the standard first Sonine estimates for the shear viscosity and the self-diffusion coefficient 
are also partially corrected by the modified approximation. Finally, the extension of the modi- 
fied first Sonine approximation to the transport coefficients of the Enskog kinetic theory is presented. 

Keywords: Granular gases; Boltzmann kinetic theory; Navier-Stokes transport coefficients; So- 
nine approximation 
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I. INTRODUCTION 

The usefulness of kinetic theory tools to describe the 
dynamical properties of granular fluids has been widely 
recognized [l|, |2|- The primary difference from normal 
fluids lies in the inelastic character of collisions, what in- 
troduces features not present in ordinary matter, such as 
the absence of equilibrium states, the spontaneous for- 
mation of clusters, and the development of high-energy 
tails, among others. The prototype model of a granular 
gas is a system composed by a large number of smooth 
hard spheres colliding inelastically with a constant coef- 
ficient of normal restitution < a < 1. In the dilute 
limit, the Boltzmann kinetic equation, suitably modified 
to incorporate inelasticity, provides a convenient frame- 
work to investigate some of the most relevant properties 
of granular gases. 

One of the main applications of the inelastic Boltz- 
mann equation is the derivation of the constitutive equa- 
tions for the stress tensor and the heat flux in a hydro- 
dynamic description. These constitutive equations dcflne 
the relevant Navier-Stokes (NS) transport coefficients 
and are determined by means of the Chapman-Enskog 
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(CE) method to solve the Boltzmann equation for suf- 
ficiently long space and time scales. While in the elas- 
tic case the velocity distribution function / is expanded 
about the local equilibrium distribution function Jm , the 
reference state for granular gases is the local version /'^"^ 
of the so-called homogeneous cooling state (HCS). The 
first-order stage of the CE expansion allows one to ex- 
press the transport coefficients in terms of the solutions 
of linear integral equations 0, II, H, la B ■ ^"^ these equa- 
tions the distribution /^^^ appears explicitly in the inho- 
mogeneous terms and implicitly through the linearized 
Boltzmann collision operator. As happens in the elastic 
case, the solutions of the integral equations are not known 
exactly, so that approximations must be introduced in 
order to get the transport coefficients as explicit non- 
linear functions of a. The standard method consists of 
approximating the solutions by the Maxwell-Boltzmann 
distribution /m times truncated Sonine polynomial ex- 
pansions. The simplest possibility is the first Sonine ap- 
proximation, where only the lowest Sonine polynomial 
is retained. The resulting expressions for the effective 
collision frequencies associated with the transport coef- 
ficients have an explicit dependence on a (due to the 
collision rules) as well as an implicit one through a linear 
dependence on the fourth cumulant 02 of f^^' . 

The reliability of the standard first Sonine approxima- 
tion has been tested in the last few years H, M: Il0> HB 
|i3,llJ,llJ,[l5, 16, 17] by comparison with computer sim- 
ulations of the Boltzmann equation by means of the di- 
rect simulation Monte Carlo (DSMC) method 1^. The 



comparisons show that the shear viscosity 77 and the self- 
diffusion coefficient D are accurately estimated by the 
first Sonine approximation, even for strong dissipation. 
The two transport coefficients k, and ^ characterizing the 
heat flux are well described by the first Sonine approxi- 
mation for moderate and small inelasticity (say a -^ 0.7). 
However, recent studies [ij; US show that the first So- 
nine approximation dramatically overestimates n and /j, 
for high dissipation {a < 0.7). 

Although the range a ^ 0.7 encompasses the region 
of practical and experimental interest, especially if one 
considers the inherent couplin g b etween gradients and 
dissipation in steady states J19l |2(]| , it is important from 
a fundamental point of view to understand the origin 
of those discrepancies for k and /i and propose alter- 
native approximations. Since in the simulations carried 
out in Refs. |14L Il5l| the transport coefficients were ob- 
tained from two-time correlation functions by means of 
Green-Kubo (GK) relations, the discrepancies might be 
due to velocity correlation effects outside the domain of 
the Boltzmann equation. However, as discussed in Ref. 
[13 , the disagreement seems to be directly related to the 
failure of the Sonine expansion truncated after the first 
term to capture the velocity dependence of the NS dis- 
tribution function. One possibility of improving the ap- 
proximation could be to consider higher order terms in 
the Sonine expansion [l^l- However, the involved algebra 
would be rather intricate and it is not obvious that the 
improvement would be significant. 

The aim of this paper is to implement an alternative 
route to the standard first Sonine approximation. The 
idea is based on the assumption that the isotropic part 
of the NS velocity distribution /^^^ is mainly governed by 
the HCS distribution /(°) rather than by the Maxwellian 
distribution Jm H- More specifically, our modified first 
Sonine approximation has the same form as the standard 
one, except that the weight function /m is replaced by 
/("^ As a consequence, the effective collision frequencies 
derived from the modified approximation have the same 
structure as those derived from the standard one, except 
that the respective coefficients of 02 differ markedly in 
both approximations. Since the high-velocity population 
in /'^"^ is larger than in /m, it is reasonably expected 
that the former captures better the influence of the high- 
velocity tail of Z*-^-' on the NS transport coefficients, es- 
pecially those related to the heat flux. In fact, the results 
show that, in the region a ^ 0.7, the modified values for 
the collision frequencies are larger than their standard 
counterparts, so that the corresponding transport coef- 
ficients are smaller in the modified first Sonine approxi- 
mation than in the standard one. As will be shown later, 
the modified estimates for the transport coefficients k 
and /i compare quite well with available computer simu- 
lations, even for extreme dissipation, in contrast to what 
happens with the standard estimates. For the remaining 
transport coefficients 77 and D, which are already well 
described by the standard approximation, the modified 
approximation provides even better values. 



The plan of the paper is as follows. The expressions of 
the NS transport coefficients obtained by the application 
of the standard first Sonine approximation are recalled 
in Sec. ^ Our modified first Sonine approximation is 
described and discussed in Sec. lIIII Next, the two Sonine 
approximations are compared in Sec IIVI with available 
and new simulation data for the transport coefficients 
K, /x, 77, and Z), both for two- and three- dimensional 
systems. The extension to the transport coefficients pro- 
vided by the Enskog theory is presented in Appendix IBl 
The paper is closed with some concluding remarks in Sec. 

m 



II. STANDARD FIRST SONINE 
APPROXIMATION 

We consider a granular gas of smooth, inelastic hard 
spheres (in d dimensions) of mass tti, diameter a, and 
coefficient of restitution a. In the low-density regime, 
the one-particle velocity distribution function /( r,v,t) 
obeys the (inelastic) Boltzmann equation |2ll |2^ . Un- 
der the conditions of weak hydrodynamic gradients, the 
CE method [23 provides a solution of the Boltzmann 
equation based on an expansion / — /(°) -f f^^^ + • • • , 
where /(°) is the local version of the HCS 21, "24]. The 
first-order distribution f^^' has the form ^, ^] 

/(^^ (V) = ^(V) • V InT + B(V) • V In n -F C,j (V) V,Uj , 

(2.1) 
where n, T, and u are the number density, granular tem- 
perature, and flow velocity, respectively, and V = v — u 
is the peculiar velocity. The functions .4,(V), B(V), 
and Cij(V) are the solutions of a set of linear integral 
equations. While >t(V) and Cijiy) obey autonomous 
equations, the equation for B(V) requires the knowl- 
edge of ^(V). However, the combination ^'(V) = 
A(V) — iB(V) satisfies a closed equation TS]- From the 
solutions to those linear integral equations, the shear vis- 
cosity 77 (associated with the pressure tensor), the ther- 
mal conductivity k, and a new transport coefficient /i 
(the two latter associated with the hcat-fiux) are formally 
given by Hi 113 



1 



v^m- 



hC 



Kq- 



d-1 l + 2a2 



2T, ,, 

n 



2C*' 



1 1 



Ko- 



|a2 



c* 



(2.2) 
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where, 770 and kq are the elastic values (in the first Sonine 
approximation) of the shear viscosity and thermal con- 
ductivity [23j, respectively [cf. Eq. ljB10|l ]. In addition, 
Q* = C(°)/i/o is the reduced cooling rate of the HCS, 



where i/q = nT/riQ is an effective collision frequency [cf. 
Eq. (|1J12|| ] , and 02 is the fourth velocity cumulant of f^^'' . 
Furthermore, in Eqs. H2.2|l - H2.4|l the (reduced) effective 
collision frequencies are 



/rfvAj(v)/:c»j(v) 

z.o/rfvA,(V)C,,(V)' 



(2.5) 



J d^rS{V) ■ CA{V) 



JdvS{V)-CA'iV) 



" 1^0 J dvS{V) ■ AiV)' " iyoJdvS{V)-A'{V)' 

(2.6) 
where C is the linearized Boltzmann collision operator 
and we have introduced the polynomials 



(2.7) 



D,j{V)^m[V,Vj--V'S,j 



S(V) 
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d + 2. 



T V. 



(2.8) 



Another relevant transport coefficient for a single gas is 
the self-diffusion coefficient D, which measures the diffu- 
sion of tagged particles in a fluid of mechanically equiva- 
lent particles in the HCS. Application of the CE method 



leads to fr = {n,/n)f(^°'> and 



/ii)(V)=X>(V).Vn, 



(2.9) 



where Ug and /s are the number density and velocity dis- 
tribution function of the tagged particles, respectively. 
As before, the function "DCV) is the solution of a linear in- 
tegral equation. The CE expression for the self-diffusion 
coefRcient D is [9l[l3| 



02 = 



16(l-a)(l-2a2) 



24d-a(57 



2(l-a)a2 



(2.13) 



However, the determination of the collision frequencies 
v*, v*^,, v* and v*j^ is much more complicated since it 
requires the knowledge of A(V), A'(V), Cy(V), and 
X>(V), respectively. These functions are the solutions of 
linear integral equations in which the HCS distribution 
f^^i appears explicitly in the inhomogeneous terms and 
also implicitly through the linearized Boltzmann opera- 
tor. From that point of view, those collision frequencies 
are Junctionals of /'°'. In order to get explicit expres- 
sions for the dependence of the transport coefficients on a 
one has to resort to some approximations. As in the elas- 
tic case, the simplest approximation consists of truncat- 
ing the Sonine polynomial expansions of A(V), ^'(V), 
Cy(V), and X'(V) after the first term. More explicitly, 
the first Sonine approximation is 



A{Y) 
A'(V) 
Cy(V) 
V{Y) 



fAliV) 



c.S(V) 

c.'S(V) 

c^Aj(V) 



where the Maxwellian 
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(2.15) 



is the weight factor in the scalar product with respect 
to which the orthogonal polynomials are defined. The 



coefficients 



and CD are directly related to the 



transport coefficients by 



d + 2nT^ 
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where Dq is the elastic value of the self-diffusion coeffi- 
cient (in the first Sonine approximation) [cf. Eq. ljBll|l ] 
and 



T] = —nT 



--rjj 



(2.18) 



'D 



/dvV.£,P(V) 
j/o/dvV-X>(V) 



(2.11) 



is the (reduced) collision frequency associated with 
the self-diffusion coefficient, Ls being the Boltzmann- 
Lorentz operator. 

So far, the expressions for the NS transport coefficients 
are formally exact, but their a-dependence through the 
quantities C*, 02, v^^ ^*, ^*,, and v*j^ is not explicitly 
known. The two first quantities (C* and 02) depend on 
the HCS distribution and are very accurately estimated 
bylHllllli 



C-^(-«^ 



l + ^«2,, 
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D 



nT 



-CD- 



(2.19) 



Note that the approximations (|2.14|) imply that A'(V) ex 
A(V), what leads to v*^, = v*^. 

Inserting the approximations H2.14|l into Eqs. (|2.5|l . 
(|2.6|) . and H2.11|l . and taking the Sonine approximation 
for /'■°-', one can evaluate explicitly the (reduced) colli- 
sion frequencies H2.5|l . H2.6|l . and H2.11|l . The results are 
iii 
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where S(V) and Dij (V) have the same polynomial struc- 
ture as S(V) and £'y(V), respectively, but must be cho- 
sen to preserve the solubility conditions |23,|23|- A simple 
calculation yields 



d + 2, , / 1 

^^(1 + "H'~32«^ 



(2.22) 



Therefore, the NS transport coefficients in the standard 
first Sonine approximation are given by Eqs. (|2.2|l - (|2.4|) 
and H2.10|l with the collision frequencies given by Eqs. 
(I2.2()(l - (|2.22(l . In addition, the coohng rate C* and the 
fourth cumulant a2 are given by Eqs. H2.12(l and H2.13|l . 
respectively. 

As said in the Introduction, while the results obtained 
in the first Sonine approximation for the shear viscosity 
11 El 111 113 and the diffusion coefficient [i|ll|l3 com- 
pare quite well with computer simulations over a wide 
range of inelasticities, the coefficients k and ^ associated 
with the heat flux show important discrepancies with 
simulation data for strong inelasticity ,1^ J^ JJJ . 



A,(V) = A,(V) 



ijV '' h 



S(V) = S(V) - ^^as^V. 



(3.2) 



(3.3) 



As before, ^'(V) ex >l(V), so that j/*, = i^* in the modi- 
fied first Sonine approximation also. The coefficients c„, 
Ck', Cfj, and cd are related to the transport coefficients 
by 



d + 2 nT'^ ( d + 8 d + 2 ., d + 4 



2 m 



-0.2 ^a2 ^^3 ] Ck, 

(3.4) 



(3.5) 



III. MODIFIED FIRST SONINE 
APPROXIMATION 



T] ^ -nT"^ (1 + a2)c^. 



(3.6) 



One of the possible sources of discrepancy between the 
standard first Sonine approximation for the transport co- 
efficients associated with the heat flux and computer sim- 
ulations could be due to the existence of non-Gaussian 
features. Although the Maxwellian distribution /m is 
a good approximation to J^^> in the region of thermal 
velocities relevant to low-order moments (hydrodynamic 
quantities), quantitative discrepancies between both dis- 
tributions are expected to be important in the case of 
higher velocity moments, such as the heat flux. The de- 
parture of j^^' from /a/ is partially accounted for by a2. 
However, in the approximation (|2.14l) the behavior of /*^^^ 
is assumed to be essentially dominated by the Maxwellian 
distribution /m- From that point of view, one might say 
that a certain mismatch exists in the standard first So- 
nine approximation applied to inelastic gases. This could 
be fixed by incorporating more terms in the Sonine poly- 
nomial expansion [13, but this would be at the expense 
of significantly increasing the technical difficulties of the 
method. 

Here we follow an alternative route, similar to the one 
discussed in Ref. |3- Specifically, we keep the structure 
of (|2.14|) . except that the distribution /'^°) is chosen in- 
stead of the simple Maxwellian form /a/ as the convenient 
weight function. According to these arguments, we take 
the approximations 
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In Eq. (|Si3J), 03 is the sixth cumulant of /(°) 24]. Note 
that in Eqs. H3.2|l - (|3.7|) no explicit form for /("5 has been 
needed to be assumed. 

In the modified first Sonine approximation, the colli- 
sion frequencies are obtained from Eqs. 12.5|l . (|2.()|) . and 
(|2.11f) by inserting the approximations l|3.1|l . and neglect- 
ing 03 and nonlinear terms in 02 . After lengthy algebra 
[Hi, one gets 
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v^ = v^., = 



d — 1 3 , , 
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A('l_a + ^dVl + a)('l + ^a2), (3.9) 
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Thus, the NS transport coefficients in the modified ffi'st 
Sonine approximation are given by Eqs. H2.2|l - H2.4|l and 
(|2.10|) . with the collision frequencies given by Eqs. (|3.8|l - 



Comparison between the standard approximations, 
Eqs. (E2ni)-(E!221), and the modified ones, Eqs. (|^ - 
(I3.1()() . shows that they only differ in the coefficient of 
the term linear in 02- In the standard approximation, 
the dependence of the collision frequencies on 02 only 
arises from the presence of the HCS distribution /'"^ in 
the linear operators £ and Cs ■ On the other hand, in the 
modified approximation there exist additional contribu- 
tions arising from the weight factor /(°) in Eq. H3.1|l and, 
in the case of v'^, also from the modified Sonine polyno- 
mial S(V). These additional contributions give rise to 
a renormalization of the coefficients of 02, which change 
dramatically with respect to their values in the standard 
approximation. More specifically, the coefficient of 02 
in Eq. (|3.8(l is at least 46 times larger than the coeffi- 
cient in Eq. H2.20|l . both being positive. In the cases of 
V* and J/^, the coefficients are negative in the standard 
first Sonine approximation, while they are positive in the 
modified Sonine approximation. Moreover, the magni- 
tudes of the coefficients in the latter approximation are 
14 and 6 times larger, respectively, than in the former 
one. These discrepancies are not significant as long as 
the magnitude of 02 is relatively small. This is what 
happens for a ^ 0.7. However, for larger inelasticity, 
the fourth cumulant 02 is not negligible |25l l2^ uS\ yOl ■ 
Since 02 > for a < 0.7, then the standard estimates for 
the collision frequencies are smaller than their modified 
counterparts. Consequently, the associated transport co- 
efficients are larger in the standard approximation than 
in the modified one. The fact that these effective collision 
frequencies associated with the transport coefficients are 
larger in the modified approximation than in the stan- 
dard one is possibly due to the overpopulation of /*•"■' 
with respect to /m for high velocities. Since the colli- 
sion rate for hard spheres increases with velocity, this 
overpopulation yields a more efficient average collisional 
transfer of momentum and energy. 

In principle, Eq. (|3.1|l can be seen as the first-order 
approximation in a polynomial expansion. For instance, 
in the case of ^(V) one can write 



^(V) 



Kd/2) 



/(")(V)V^c,Lr^(c^), 



(3.11) 



k=l 



where {Lj, (x)} is a set of orthogonal polynomials with 
respect to an inner product involving f^^'. If f^^' is re- 
placed by Jm, then the polynomials Lj, {x) become 
the generalized Laguerre polynomials L^, [x) . The ex- 
pansion (|3.11() differs from the one considered in Ref. Q 

in the use of the modified polynomials Lp. (x) instead 
of the conventional polynomials Lj^ {x)i which do not 
constitute an orthogonal set in this case. The recursive 

procedure to get the polynomials Lf, (x) in terms of 
the cumulants of f^^' is briefly described in Appendix 

HI 



IV. COMPARISON WITH COMPUTER 
SIMULATIONS 



In this Section we compare the theoretical expressions 
for the transport coefficients obtained from the standard 
and modified first Sonine approximations with available 
and new computer simulations. 



A. Heat flux 

The NS transport coefficients associated with the heat 
fiux are k and /i. These transport coefficients have 
been measured from the GK relations l31| bv means of 
the DSMC method fl§\, both for two- gj and three- 
dimensional jlSj systems. In addition, the coefficient 
k' = K — n/i/2T has been measured in DSMC simula- 
tions by an alternative method based on the application 
of an external force in the three-dimensional case [13I ■ 

Figures Hand 121 show the a-dependence of the reduced 
transport coefficients k/kq and fj,/(TKQ/n), respectively. 
It is apparent that the standard first Sonine approxi- 
mation significantly overestimates both transport coef- 
ficients for strong inelasticity. On the other hand, the 
modified approximation compares well with computer 
simulations, even for low values of a, especially in the 
three-dimensional case. This refiects the fact that the 
modified approximation is more accurate than the stan- 
dard one in describing the effective collision frequencies 
for heat transport. 

Since both k and /i are overestimated by the standard 
approximation, it could happen that, by a cancelation 
of errors, the transport coefficient k' — k — nfi/2T might 
be well captured by that approximation. However, this is 
not the case. The comparison between the computer sim- 
ulations for k' obtained from the two alternative methods 
of Refs. [13 and [l7| and the two theoretical approaches 
is shown in Fig. |5l As in the cases of k and /i, the mod- 
ified Sonine approximation for k' agrees quite well with 
the simulation results. 

The good agreement found in Figs. HEl between the 
modified first Sonine approximation and the simulation 
data for the heat flux transport coefficients suggests that 
the NS distribution functions -4.(V) and A-'CV) are well 
represented by the forms H3.ll) . To test this expectation, 
we compare now the standard and modified Sonine ap- 
proximations for A. (V) with simulation data presented 
in Ref. [I3I for the three-dimensional case. By symmetry 
arguments, the function jA. (V) can be written as 



A!(V) = Xv^'fM{VMc')y, 



(4.1) 



where A ~ l/v27rncr^ is the mean free path, vq — 
y/2T/rn is the thermal speed, and <I>(c^) is a dimension- 
less isotropic function of the scaled velocity c = V/vq. 
All the information contained in $(c^) is retained by the 





FIG. 1: (Color online) Plot of the (reduced) thermal conduc- 
tivity k/ko as a function of a for hard disks (top panel) and 
hard spheres (bottom panel). The dashed and solid lines rep- 
resent the standard and modified first Sonine approximations, 
respectively. The symbols are DSMC results obtained from 
the GK relations [TllTsl. 



FIG. 2: (Color online) Plot of the (reduced) transport coeffi- 
cient fj,/{KoT/n) as a function of a for hard disks (top panel) 
and hard spheres (bottom panel). The dashed and solid lines 
represent the standard and modified first Sonine approxima- 
tions, respectively. The symbols are DSMC results obtained 
from the GK relations |l4 l ll5 || . 



marginal distribution IT3 
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According to the standard approximation H2.14|l . 
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In contrast, the modified approximation (|3.1|) yields 



(4.4) 



*(c^ 



4k' 
5nXvo 

4k' 
5nXvo 



1 



a2 



5c' 



15 



(l + a2)-c2 



1 + ^02 



c + a2 , 

2 V 16 



105 21 



'-^c^ + lc^ 



(4.5) 



where in the last equality we have neglected nonlinear 
terms in 02 . The corresponding marginal distribution is 
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FIG. 3: (Color online) Plot of the (reduced) modified thermal 
conductivity k'/kq as a function of a for hard spheres. The 
dashed and solid lines represent the standard and modified 
first Sonine approximations, respectively. The symbols are 
DSMC results obtained from the GK relations IirII (triangles) 
and from the application of an external force Jl7| (circles) . 



The function (p{c'^) is plotted in Fig.^jfor a = 0.5 and 
a = 0.3. We see that the modified first Sonine distribu- 
tion (|4.6|) captures reasonably well the main features of 
the true distribution, especially for a = 0.5. On the other 
hand, the standard first Sonine distribution 14.41) strongly 
disagrees with the simulation data for c^ > 6. This ve- 
locity region has a significant infiuence on the evaluation 
of the thermal conductivity at high inelasticity [I3I . 




velocity correlations in the correlation function involved 
in the GK expression of the shear viscosity. These veloc- 
ity correlations are larger than in the case of the heat-flux 
transport coefficients |l4J . 
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FIG. 4: (Color online) Plot of the marginal distribution func- 
tion (p{c^) obtained from DSMC (solid lines) for a — 0.5 (top 
panel) and a — 0.3 (bottom panel). The dashed and dotted- 
dashed lines represent the standard and modified first Sonine 
approximations, respectively. 



B. Pressure tensor 

Although the a-dependence of the shear viscosity is 
well described by the standard first Sonine approxima- 
tion, it is worthwhile comparing the modified first So- 
nine estimate against computer simulations. To the best 
of our knowledge, the NS shear viscosity has been mea- 
sured in DSMC simulations by three alternative methods: 
(i) by analyzing the time decay of a weak transverse shear 
wave in the HCS Jj; (ii) from the GK relation [11 [TH: 
and (iii) by the application of a homogeneous external 
force |IQ|. The simulation data obtained by these meth- 
ods and the two theoretical approximations are presented 
in Fig.[S| The data obtained from the method (i) are re- 
stricted to d = 3 and a > 0.7 [a|, while the ones from 
the method (ii) are available for d = 2 [Ij] and d = 3 
|15| . Regarding the method (iii), the data correspond- 
ing to d = 3 for a > 0.6 were reported in Ref. [iq . 
while those corresponding to d = 3 for a > 0.5 and to 
d = 2 have been obtained in this work. We observe that 
up to Of ~ 0.6 the simulation data are consistent among 
themselves and also with both theories. However, for 
higher inelasticities, there is a certain discrepancy (less 
than 10%) between the data reported in Refs. [H Il5j 
and those presented here, the former being close to the 
standard estimates and the latter being close to the mod- 
ified estimates, especially in the three-dimensional case. 
The small difference between our simulation data and 
those of Refs. [ijj UM might be due to the influence of 
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FIG. 5: (Color online) Plot of the (reduced) shear viscosity 
rj/rjo as a function of a for hard disks (top panel) and hard 
spheres (bottom panel). The dashed and solid lines repre- 
sent the standard and modified first Sonine approximations, 
respectively. The symbols are DSMC results obtained from 
the decay of a sinusoidal perturbation [g (squares), from the 
GK relations |14lll5l (triangles'), and from the application of 
an external force \W^ (circles). In the latter case, the data 
corresponding to d = 3 for a < 0.5 and to d = 2 have been 
obtained in the present work. 



In conclusion, while the standard first Sonine approxi- 
mation does quite good a job for the shear viscosity, it is 
fair to say that the modified approximation is still better, 
especially for three-dimensional systems. 



C. Self-diffusion 

Finally, we consider the self-diffusion coefficient. This 
coefficient has been measured in computer simulations 
from the mean square displacement of a tagged particle 
in the HCS [3, [ifl UM i ^^ ^^^^ ^^ tiy the decay of a sinu- 
soidal perturbation in the concentration of tagged parti- 
cles 9] . We observe in Fig. that both Sonine approxi- 
mations provide a general good agreement with simula- 
tion data. However, the standard approximation slightly 
overestimates the self-diffusion coefficient at high inelas- 
ticity, this effect being corrected again by the modified 
approximation. 
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FIG. 6: (Color online) Plot of the (reduced) self-diffusion coef- 
ficient D/Dq as a function of a for hard disks (top panel) and 
hard spheres (bottom panel). The dashed and solid lines rep- 
resent the standard and modified first Sonine approximations, 
respectively. The symbols are computer simulation results ob- 
tained from the mean square displacement. In the top panel, 
the triangles correspond to molecular dynamics results |9| and 
the circles are DSMC results obtained in the present work. In 
the bottom panel, the triangles are DSMC results reported in 
Ref. [31 and the circles are DSMC results presented in Ref. 
[13 {ct > 0.6) and obtained in the present work (a < 0.5). 



V. CONCLUDING REMARKS 

This work has been mainly motivated by the disagree- 
ment found at high dissipation between the simulation 
data for the heat flux transport coefficients and the ex- 
pressions derived from the standard first Sonine approxi- 
mation fl4l Il5lll7l | . Although this disagreement appears 
beyond the region of inelasticity of practical interest, it is 
physically relevant from a fundamental point of view to 
propose alternative theoretical approaches that correct 
the limitations of the standard approximation. Here, we 
have implemented a modified version H3.1|l of the first 
Sonine approximation (|2.14() . where the weight function 
is no longer the Maxwell-Boltzmann distribution /m but 
the HCS distribution f'^^\ Moreover, in order to preserve 
the solubility conditions, the polynomial S(V) defined by 
Eq. H2.8|l must be replaced by the modified polynomial 
S(V) defined by Eq. |j23- The idea behind the modified 
method is that the deviation of /'■*'•' from /j\/ has an im- 
portant influence on the NS distribution f'^^\ so that the 
latter is better captured by the approximation H3.1|l than 
by the approximation (|2.14|) . In other words, the rate 
of convergence of the polynomial expansion is expected 
to be accelerated when /(°) rather than Jm is used as 



weight function. 

The structure of the transport coefficients is common 
in both approximations. They are given by Eqs. (|2.2|l - 
(|2.4|) and H2.10|l . However, the a-dependence of the char- 
acteristic collision frequencies differs in both methods. In 
the standard first Sonine approximation, those collision 
frequencies are given by Eqs. (|2.2U|I - 12.22|I . while they are 
given by Eqs. (|3.8|I - H3.1UII in the modified first Sonine ap- 
proximation. It is apparent that the distinction between 
both approximations occurs in the value of the coeffi- 
cient of the fourth cumulant 02. In the standard approx- 
imation that coefficient conies from the dependence of 
the linearized Boltzmann collision operator on f^^' only, 
while in the modified approximation it also comes from 
the assumed form for f^^\ The effect of the latter contri- 
bution becomes more important than that of the former, 
so that, for each collision frequency, the coefficient of 02 
changes dramatically from the standard approximation 
to the modified one. 

As observed in Figs. EOI the modified approximation 
significantly improves the a-dependence of k, /i, and the 
difference k' — k — nfj,/2T. This is the primary result of 
this paper. Additionally, as shown in Figs. [S] and EJ the 
slight discrepancies between simulation and the standard 
first Sonine estimates for the shear viscosity j] and the 
self-diffusion coefficient D are partially corrected by the 
modified approximation. 

Although the results reported here have been restricted 
to a low-density granular gas described by the inelas- 
tic Boltzmann equation, they can be straightforwardly 
extended to finite density in the framework of the En- 
skog kinetic theory. In that case, application of the CE 
method shows that the effective collision frequencies are 
the same as in the dilute limit, except for a density- 
dependent factor X 13, Q- The explicit expressions for 
the NS transport coefficients are presented in Appendix 
IbI We expect that these Enskog results can stimulate 
the performance of molecular dynamics simulations to 
test whether or not the modified first Sonine approxi- 
mation improves again over the predictions of the stan- 
dard approximation, especially in the case of the heat 
flux transport coefficients at high inelasticity. 
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APPENDIX A: MODIFIED POLYNOMIAL 
EXPANSION 



coefficient as in Lj^' (x) . The orthogonalization condition 
(4^\l[^^) = for ^ < fc - 1 gives 



Given a weight function 



w{x) 



1 + Y.^uL^' \^) 



k=2 



(Al) 



the mathematical problem consists of finding a set of 

polynomials {Lj, (a;)} such that they are mutually or- 
thogonal with respect to the scalar product 



(/i,/2)-/ dxxPw{x)h{x)f2{x). (A2) 

If flfe = 0, then w{x) = e^^ and one has the Laguerre 

polynomials, i.e., Lf. {x) = Lj^ (x). In the general case, 

the polynomials Lf. (x) can be obtained following the 
Grani-Schniidt orthogonalization procedure. Suppose 

the polynomials L^ {x) with £ < fc— 1 are already known. 
The next unknown polynomial can be written as 



^k W — ^k ^ 



fc-1 



+ VcfLr(x), 



i=0 



(A3) 



Jfc) 



(fc) [Lj ,x'') 



-rip), 



(A4) 



This closes the construction of Lj, (x) and the process 

can be recursively continued. Since Li^_^{x)x''~^^ is a 
polynomial of degree 2k, and given the orthogonality 
properties of the Laguerre polynomials L^~ (x) appear- 
ing in the representation (|A1|I . it is straightforward to 
see that only the first k coefficients {a£,£ < k} appear 

in Lj. (x). The norm of Lf. (x), however, involves the 
coefficient a2k+i- 



The three first polynomials are Lq (x) — 1 



-(p) 



where 



LY'{x) = {p+l){l + a2)-. 



Kp) 



(A5) 



ix)^'-x^ + cl^^+cr^tfix), (A6) 



where the coefficients c[ are to be determined. One of 
them can be fixed by the standardization condition. For 
instance, we can take cj^ = (— l)''"/fc!, which is the same 



(2) ^_±^P±lt,^o ^ /47^ 

"o = 2 (l + Saa-ag), (A7) 



^(2) _ p + 22 + 2ip + 7)a2 - (p + l)a2(3Q2 - as) - (3p + ll)a3 + jp + 3)a4 
^^ ^ 2 1 -^ (p -^ 4)a2 - (p -h l)a2 - (p + 2)a3 

I 



(A8) 



APPENDIX B: TRANSPORT COEFFICIENTS 
FOR A DENSE GRANULAR GAS 

In this Appendix, we give the expressions for the NS 
transport coefficients obtained from the Enskogkinetic 
equation by the application of the CE method |j, Q in 
the first Sonine approximation. 

The bulk viscosity (which vanishes in the dilute limit) 
is 



7 = '70 



where 



22rf+l / ^ 



,rf/2 



2^^(1 + ^/2) 



(Bl) 



(B2) 



is the solid volume fraction and x((/>) is the pair corre- 
lation function at contact. The shear viscosity rj has a 



kinetic part rj'^ and a coUisional part rj'^, where 



m 

X 



< - 2^' 



■)d-2 



1- 



d + 2 



(l + a)(l-3a)0x 



Td-l 



d + 2 



0X(1 + «)»?'= 



d + 2 



(B3) 



(B4) 



Analogously, the coefficients associated with the heat 
flux have also kinetic and coUisional contributions. They 
are 



Ko d — \ 
X d 

2d-3 

+3 



d+2 



iiy:-2cr'[i + 2a2 

0X(1 + af [(1 + a)a2 - 1 + 2a]|, 



(B5) 
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3^0x(l + «)-* 



22d+l(^_l) 



{d + 2)^Tr 



M 



X I 1 + j^a2 I Ko, 



T Ko 



rx(i + «) 



(B6) 



n X 



Ho 



<c 



-a2 



-iia^ - 3a + 10 + 2d)a2 - ail - a) 
6 



id-2 



Ai" = A* 3 



<i+2 



0x(l + a)- 



(B7) 



(B8) 



In Eq. IJB7|) . we have introduced the quantity ^ = 
FinaUy, the self-difFusion coefficient is simply given by 



D = 



Dad + 2 1 



L A* 



(B9) 



In the above equations, t^o, «^0j and Dq are the elastic 
values of the low-density shear viscosity, thermal conduc- 



tivity, and self-diffusion coefficient, respectively. They 
are given by 



_nT _ d{d + 2) 770 

vq Z[d — 1) m 



Dn = 



2d T 
d + 2 rnvQ ' 



(BIO) 



(Bll) 



where the collision frequency vq is defined by 



i^O 



8 T:(d-i)/2 
d + 2 r(d/2) 



■ncT 



d-l 



1/2 



(B12) 



In addition, the reduced cooling rate C* for a dilute gas 
in the HCS is given by Eq. (|TT2|) . 

The expressions (|B3|) ~ (IB9|I are common to the stan- 
dard and the modified first Sonine approximations. How- 
ever, as discussed in the main text, they differ in the a- 
dependence of the (reduced) effective collision frequencies 
V*, V*, and i/|), which are given by Eqs. H2.20|l - H2.22|l in 
the standard approximation and by Eqs. (|3.8f) ~ l|3.10|) in 
the modified one. Note that, since the coUisional contri- 
butions depend on their kinetic counterparts, they also 
differ in both approximations. 
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